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NATIONAL ADVISORY COMMITTEE FOR AERONAUTIC 
TECHNICAL NOTE 4207 

EFFECT OF A STRINGER ON THE STRESS CONCENTRATION 
DUE TO A CRACK IN A THIN SHEET 



By J. Lyell Sanders , Jr. 


SIMMARY 


A coefficient is obtained for determining the effect of a reinforcing 
stringer on the stress concentration factor at the tip of a crack in a 
thin sheet. The results are given for the case in which the stringer Is 
intact and for the case in which the stringer is broken. In the first case 
the stress concentration factor -for the stringer is also given. 


INTRODUCTION 


Seme damage to aircraft structures due to fatigue or accident is sta- 
tistically inevitable ; thus, the fail-safe concept has entered into design 
considerations. One of the problems associated with this concept is the 
determination of the static strength of cracked parts. The mechanism of 
static failure of a structure weakened by the presence of a crack Is by 
no means completely understood at the present time. However, an engi- 
neering theory which seems to hold some promise has recently become 
available (ref. l) . In this theory the significant quantity determining 
the strength of the cracked structure is the stress concentration factor 
at the end of the crack (corrected for plasticity and the so-called size 
effect) . The fundamental Information needed to apply the method Is the 
stress concentration factor obtained from elasticity theory. 

For many configurations an exact solution for the stress distribu- 
tion from the theory of elasticity is very difficult to obtain. However, 
a considerable amount of information that is useful and adequate for 
practical applications has been obtained by making various idealizations 
and simplifications of the problems . As a further contribution, the 
results contained in the present paper were obtained. 

The problem considered in the present paper is the determination of 
the relieving effect of a reinforcing stringer on the stress concentration 
at the tip of a crack in a thin sheet. The crack runs perpendicular to 
the stringer and extends an equal distance on either side of It. The 
state of stress in the sheet far away from the crack is a tensile stress 
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parallel to the stringer. The stress concentration factor for a crack in 
a thin sheet may he determined from a known formula. (See ref. 1.) The 
factor by which this known result can be multiplied in order to correct 
for the presence of the reinforcing stringer is determined in the present 
paper. The stress concentration factor in the stringer due to the crack 
and the correction factor for the crack in the case in which the stringer 
is broken are also found -in the analysis. The results are presented 
graphically and in tabular form. 

SYMBOLS 

A 
B 
b 

C ratio between stress concentration factors for a cracked sheet 

with and without a stringer 

C* 

E 

E st 

F analytic function defined in equation (ll) 

G shear modulus for sheet material 

H stress function (see eqs. (2)) 

Iq, I]_ modified Bessel functions of first kind 

Kq, K-]_ modified Bessel functions of second kind 

Lq, L-j_ Struve functions of imaginary argument 

P load in stringer at its intersection with crack 

F load concentration factor for stringer, PE/cAE st 

R( ), l( ) real part of and imaginary part of 


corresponds to C in case where stringer is broken 
Young's modulus for sheet material 
Young's modulus for stringer material 


cross-sectional area of stringer 
function defined in equation ( 30 ) 
length of crack on one side of stringer 
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s dummy variable of integration 

t sheet thickness 

u displacement in X]_ -direction 

x, y dimensionless coordinates (see eqs. (4)) 
x l> yi physical coordinates 

z complex variable, x + iy 

7 Euler's constant, 0.57722 

£, tj complex variables 

0 dummy variable 

A similarity parameter, 2btE/AE s ^. 

cr direct stress in sheet at infinity 

cr x , r direct and shear stresses in sheet 

cr x ^ 0 direct stress in sheet with a crack but without a stringer 

$ analytic function, 0 + Ljr 

G > 0 corresponds to $ for a sheet without a stringer 

0 dimensionless stress function (see eqs. (4)) 

i|r dimensionless displacement (see eqs. (4)) 

Primes indicate differentiation and the notation ~ indicates an 
asymptotic relationship. 

ANALYSIS 


Two simplifications of the problem are made in the present analysis. 
One simplification is that the sheet is assumed inextensional in the 
direction parallel to the crack. This orthotropic sheet was introduced 
by Hildebrand (ref. 2) and greatly simplifies the equations of plane 
stress. The other simplification is to treat the crack as a straight- 
line segment and assume that the strength of the stress singularity at 
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the end of the idealized crack is a measure of the stress concentration 
due to a thin crack with a small, but nonzero, radius of curvature at its 
end. The effect of the stringer on the strength of the stress singularity 
is found by solving the two similar problems of the cracked sheet with and 
without the stringer. The desired correction factor previously defined is 
taken to be the ratio of the two strengths thus found. 


Formulation of Boundary -Value Problem 

The thin sheet with a crack and attached stringer is represented in 
figure 1. According to the orthotropic plane -stress theory of reference 2, 
the stress -displacement relations are 





*1 J 

where a x and r are the direct and shear stresses, respectively, E 

is Young's modulus, G is the shear modulus, and u is the displacement 
in the xq -direction. The displacement in the yq-direction is zero from 

symmetry. Equilibrium is satisfied if the stresses are given in terms of 
a stress function H as follows: 





( 2 ) 


An equation of equilibrium for the stringer may be obtained by con- 
sidering the portion of the stringer from the origin to xq as a free 
body. (See fig. 2.) The displacements in the stringer must be the same 
as those in the sheet along the xq-axis; hence, from equation (l) the 

E s t 

stress in the stringer must be a x , where E s ^. is Young's modulus 

E' 

for the stringer material. The required equilibrium equation is thus 
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P = avA +2 I 1 tT dx! 
x E 


r 


= A Vb 5iH_ _ f* 1 &L 


E dy n 


Jo ^ 


= A 2tH 

E ay L 


(3) 


where P is the load in the stringer at x^ = 0, A is the cross- 

sectional area of the stringer, t is the thickness of the sheet, and 
H(0,0) has arbitrarily been chosen as zero. 

Introduce dimensionless variables and parameters as follows: 


x-j = /Sbx 

1 Jo 


y x = ^ 

H = ab0 


u =-2b-t 


/EG 



(4) 


where a is the direct stress in the sheet at infinity. 

The following equations may now be obtained by eliminating ct x and 
t from equations (l) and (2): 


ci^ _ &jr_ _ cty 

Sx Sy Sy Sx 


(5) 


which are the Cauchy -Riemann equations . It follows that 
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0 + it = 3>(x + iy) = $(z) 


( 6 ) 


where $ is an analytic function of the complex variable z. 

Because of sy mme try, a boundary -value problem may be formulated 
for the upper half -plane only. From equation ( 3 ) 

- yf> = P (y = 0) (7) 

oy 


Since a x = 0 on the crack 


0=0 (x = 0, 0 ^ y $ l) (8) 


From symmetry, u = 0 on the line (x = 0, y = 1) and at the point 
(x = 0, y = 0); therefore, 

♦ - 0 - °> y - ?:>> ( 9 ) 

(x = 0, y = 0) 


Since 


and r — ^0 at infinity. 


$ ~ -iz 


(z — >~) (10) 


Solution of Boundary -Value Problem 

The boundary -value problem just formulated for $> is of the mixed 
type. However, the problem may be reduced to the more familiar Dirichlet 
type by introducing a new unknown function F defined as follows: 


F(z) = 4* + iA4> + iF 

- (I - *) ' *(1 - - F ) (11) 


The boundary -value problem for 


F is as follows: 
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II 

O 

(x = 0, y > 1) 

(12) 

1(F) = 0 

(y = 0 ) 

» 

(15) 

1(F) = F 

(x = 0, y < l)_ 



F ~ Az 

( z — >°°) 


(14) 


The solution for F is found to be 


2P , Jz 2 + 1 - 1 


F = — log 

it 


A^z 2 + 1 + 




(15) 


+ l 


where the arbitrary constant C is real and the radical is positive on 
the positive real axis. From equation (ll ) , it follows that 


+ iA<& = F - iP 


( 16 ) 


Solving for $ with the use of the relationship 0 ~ -iz at infinity 
gives 






Z 2 + 1 + 


-iAz 





(17) 


Interpretation of C . - The stress concentration factor P for the 
stringer and the coefficient C are as yet unknown. Before proceeding 
to determine them as functions of A, it is convenient to show that C 
itself is the ratio between the stress concentration factors for a 
cracked sheet with and without a stringer. 

The solution to the problem of determining the complex stress function 
for the cracked sheet without a stringer may be obtained from equations ( 15 ) 
and ( 16 ) by letting A — which is equivalent to letting E S ^A — >0 

since A = ^ The result is 
S st^ 
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$ o = -i \jz^ + 1 (l8) 

The stress field is determined from the derivative of the stress 
function, which in the neighborhood of z = i (the tip of the crack) 
behaves as follows : 

<t ' ~ (z — >i) (19) 

0 J 2i(z - i) 


There is evidently a singularity at z = i. For the cracked sheet with 
a stringer. 


~ - (z — ^i) (20) 

^2i(z - i) 

as is evident from equation (l6) since $(i) =0 from equations (8) 
and (9). On the line (x = 0, y > l), r = 0 and cr x = icr<&'. Thus, at 
the tip of the crack, 


x.o 


$ ' 
o 


~ C 


( 21 ) 


and C is evidently the required ratio. 

Determination of C and P .- ®ie two conditions available for 
determining the two unknowns P and C as functions of A are 

$>( 0 ) = 0 
3>(i) = 0 

which follow from equations (8) and (9). When applied to equation (17), 
these conditions yield 

m = -E - 1 + f°k i og El L lI + = 0 

A Ji~ + y 



(25) 
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®(i) = -E + e A 

A 



k 2 + i - 1 

t 


+ 


ILt-sV^at 


= 0 


(24) 


By use of equation (24), equation (23) becomes 

-* + ^ - i ) + r (? ios + ^ e " Ed? = ° <25) 

When the log term is integrated by parts, equations (25) and (24) 
become, respectively. 


2P 

f 0 d£ + 

r° it + c „» EdC . 0 

( 26 ) 

ijtA t 

1 iji 2 + 1 L 

4 A 2 + 1 


2P 

p 1 - ^ « ♦ 

P .« + c P«ta 5 - 0 

( 27 ) 

irtA ^ 

k ?A 2 + 1 

>4 A 2 + I 



Next make the substitution £ = Itj to obtain 

3 ? r 1 l - e~ AT| an - f 1 c ~ t >, e -^idn = 1 (28) 

nJTT"? Jo 


E .af 1 a t, + “ , c - 1 s -Hm . 0 (29) 

* **01 ^TI Ui ^2TT 


The definite integrals occurring in equations (28) and ( 29 ) are 
expressible in terms of known functions as follows: 




Kq(s ) ds 


V 
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= Kq(A) 



= k x (a) 


where Kq and K-j_ are modified Bessel functions of the second kind. 
Let 


B(A) 



(30) 


The function B(A) and its first two derivatives are expressible in terms 
of the Struve functions of imaginary argument Lq and L]_ and modified 

Bessel functions of the first kind Iq and I]_ as follows (see ref. 3) • 


B(A) 




[lo(-) " Io(s)]&3 


■\ 


B'(A) 


b"(a) 









(31) 


The first integrals of the modified Bessel functions can be expressed as 
follows : 



PA 

Kq(s)cLs 

uo 


AI 0 (A) + ^[^ 0 (A) L l(A) - Ii(A)Lq(A)] 

> 


>Kq(A) + & K 0 (A)L 1 (A) + K 1 (A)L 0 (A) 


( 32 ) 
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Equations (28) and (29) may now be written 


— B - B'C = 1 + B" 


^ I KQ(s)ds + K 0 C = K x 




( 33 ) 

( 34 ) 


where, in equation (3*0, use has been made of the formula I Kq(b)cIs = — . 

UO 2 

Equations (33) and (34) may be solved for C and P to give the following 
results 


EKq - (l + B 


C = 


") K 0 (s)ds 


BK 0 + B'J Ko(s)ds 


( 35 ) 


kA B'K l + (1+ b”)k 0 

2 Oco 

BK 0 + B* J Kq(s ) ds 


( 36 ) 


Solution for broken stringer .- For the case in which the stringer is 
broken at x = 0, the factor C must be replaced by C* obtained from 
the solution to the boundary-value problem appropriate for the broken 
stringer. The solution to this problem is easily obtained from the one 
already given. It is only necessary to set P = 0 and drop the require- 
ment <D(0) = 0* The requirement _$(i) = 0 is retained and leads to 
equation (34) as before, except P = 0. Thus, 


C* 


^1 


( 37 ) 


Computation of results . - Tables of the modified Bessel functions are 
readily available. Values of the Struve functions Lq and may be 

obtained from tables given in reference 4 for the range 0 S A = 10 at 


intervals of 0.1. No tables seemed to be available for 



Lo(s)ds; 
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therefore, values of this function for values of A ^ 2 were computed 
from a power series. For values of A between 2 and 6 , it was more con- 
venient to compute the function B by numerical integration from the 
formula 


BOO 


iV2, 


uo 


P -A sin 0) ae 


(38) 


which was obtained from equation (30) by an obvious substitution. For 
values of A > 6 , the function B was computed from the asymptotic 
series 


B(A) 


»v 


7 + log 2A - 2 



(39) 


which also was obtained from equation (30) by well-known methods 
(7 = 0.57722 is Euler's constant). For the same values of A, the 
function B' and B" were computed from the derivatives of equation (39)- 

Results of the computations for C, P, and C* for A 5 100 are 
given in table I within slide-rule accuracy. The results are also plotted 
in figures 3 to 5. _For values of A > 100 , the following asymptotic 
formulas give C, P, and C* within slide -rule accuracy: 


7 + log 2A 


p „ rt A + O.875 
27 + log 2A 

c* ~ 1 + JL 

2A 


(40) 


NUMERICAL EXAMPLE 


Consider a sheet 0.1 inch thick reinforced by a stringer made of the 
same material with an area of 0.5 square inch. A crack 6 inches long 
extends 3 Inches on either side of the stringer. The effective radius of 
curvature p e at the tip of the crack is taken to be 0.002 inch. 

According to a well-known formula, the theoretical stress concentration 
factor Krp at the tip of the crack is given by 
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% - 1 + 2j£ - 1 + 2 


2 — = 78.6 


0.002 


In the present example, 


, 2bt 2 X 3 X 0.1 _ _ 

A — — - — — — J- • <— 

A 0.5 


From table I, 


C = 0.688 
P = 1-593 


The corrected stress concentration factor K^,’ at the tip of the crack 
is thus 


IQp' = CKij * 0.688 X 78.6 = 54.1 


The stress concentration factor In the stringer is P = 1.593* The 
Neuber stress concentration factor for the crack, taking size effect 

into account, is (see ref. l) 

%*!( 1 + v ) - 2 T-5 

In practical applications, of course, this large stress concentration 
factor is considerably reduced when corrected for the effect of 
plasticity. (See ref. 1 for details.) 

DISCUSSION 


Examination of figures 3 to 5 reveals at least two qualitative 
features of the results which are of practical interest. One is the 
appreciable stress concentration in the stringer and the other is the 
detrimental influence of the stringer once it has broken. These results 
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confirm intuition. The stringer is expected to carry part of the 
load refused by the sheet because of the crack. If then the stringer 
breaks, the two intact halves of the stringer carry load into the region 
of the sheet around the middle of the crack which tends to spread the 
crack more than if there were no stringer. 

Because of the idealizations made in obtaining the theoretical solu- 
tion, some caution should be observed in applying the results. In the 
analysis the stringer is assumed to be continuously attached to the sheet 
along a line . In reality the stringer has some finite width and may be 
attached to the sheet by means of rivets . Thus the theoretical results 
cannot be expected to be accurate for crack lengths shorter than two or 
three times the rivet spacing, or two or three times the width of an 
integral stiffener. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 23 , 1957- 
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TABLE I.- COMPUTED VALUES OF C, F, AHD C* 


c 

P 

C* 

A 

C 

P 

C* 

0.637 

1.000 

OO 

1.7 

0.699 

1.806 

1.265 

.645 

1.060 

4.06 

1.8 

.701 

1.848 

1.251 

.652 

1.115 

2.73 

1.9 

.703 

1.889 

1.239 

.658 

1.168 

2.23 

2.0 

.704 

1.930 

1.228 

.662 

1.219 

1.960 

3.0 

.718 

2.32 

1.156 

.667 

1.269 

1.791 

4.0 

.729 

2.69 

1.119 

.670 

1.318 

I.676 

5.0 

• 737 

3-05 

1.096 

.674 

1.366 

1.590 

6.0 

.744 

3.40 

1.080 

.677 

1.413 

1.524 

7-0 

• 750 

3.73 

1.069 

.680 

1.459 

1.472 

8.0 

.755 

4.06 

1.061 

.683 

1.504 

1.430 

9.0 

• 759 

4.38 

1.054 

.686 

1.549 

1.394 

10.0 

.763 

4.70 

1.049 

.688 

1.593 

1.365 

15.0 

• 777 

6.21 

1.033 

.691 

I.636 

1.339 

20.0 

.787 

7.64 

1.025 

.695 

. 1.679 

1.317 

30.0 

.800 

10.35 

1.017 

.695 

1.722 

1.297 

50.0 

.816 

15.40 

1.010 

.697 

1.764 

1.280 

100.0 

.834 

27.0 

1.005 
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0.1 i v 10 100 


Figure 3 . - Variation with A of the ratio C of the stress concentra- 
tion factors in a cracked sheet with and without a stringer for the 

p*Vs4-TP 

case in which the stringer is intact. A = 

AE s t 




Figure 5.- Variation with A of the ratio C* of the stress concentra 
tion factors in a cracked sheet with and without a stringer for the 

case in which the stringer iB broken. Tv = 

A^st 
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